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Abstract 

Vehicles  moving  on  unpaved  surfaces  can  be  significant  sources  of  fugitive 
dust  emissions.  The  performance  of  vehicle  operational  characteristics  as  well  as 
evaluating  vehicle  dust  mitigation  components  including  filters,  radiators,  seals, 
and  windshields  require  models  for  the  behavior  and  transport  of  emitted  fugitive 
dust.  Much  of  the  work  on  modeling  fugitive  dust  emissions  from  vehicles  have 
focused  on  large  scale  atmospheric  transport  for  environmental  quality  purposes. 

In  this  study,  however,  we  present  a  dust  emission  and  transport  model  suitable 
for  modeling  near- field  effects  of  fugitive  dust  emissions.  The  presented  model  is 
validated  numerically  for  the  prediction  of  dust  plumes  in  a  region  near  vehicles 
through  a  careful  comparison  to  the  available  experimental  data.  This  model  pro¬ 
vides,  for  the  first  time,  a  validated  dust  emission  and  transport  models  suitable 
for  the  near  field  of  moving  vehicles.  The  weakness  of  the  current  near-field  dust 
emission  and  tranport  technology  is  also  discussed. 

key  word:  Fugitive  Dust  Emission,  Computational  Fluid  Dynamics,  Multiphase 
Flow 


1  Introduction 

Particles  suspended  in  air  by  vehicular  movement  on  paved  and  unpaved  roads  are  a 
major  contributor  to  fugitive  dust  emissions,  which  in  turn  can  have  significant  impact 
on  many  aspects  of  vehicle  performance.  Scattering  caused  by  suspended  dust  particles 
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and  dust  deposition  on  wind  shields  and  mirrors  may  can  impair  visibility,  especially  in 
convoy  types  of  settings.  Operationally,  dust  can  be  ingested  by  the  engine  and  power 
train,  increasing  maintenance  expense  and  reducing  vehicle  life-time.  In  commercial 
applications,  the  transport  of  particulates  around  a  vehicle  is  of  particular  interest 
in  determining  vehicle  soiling.  Therefore,  accurate  prediction  of  dust  emission  and 
transport  in  the  near  field  of  moving  vehicles  is  needed  to  evaluate  the  impact  of  dusty 
environments  on  many  aspects  of  vehicle  performance. 

There  are  a  number  of  challenges  involved  in  the  modeling  of  dust  emission  and 
transport  caused  by  moving  vehicles  traveling  on  unpaved  roads.  First,  the  dust  gener¬ 
ation  mechanism  is  dependent  on  complex  phenomena  that  occur  at  the  contact  between 
the  vehicle  tires  or  treads  and  the  dusty  surface.  These  phenomena  includes  the  process 
of  crushing  large  particles  into  smaller  pieces  by  the  vertical  pressure  of  the  wheel  and 
then  the  transport  of  the  crushed  particulates  by  the  tire  due  to  horizontal  friction  and, 
finally,  the  ejection  of  particulate  matter  into  the  air  due  to  centrifugal  force.  Saltation, 
which  involves  the  effects  of  large  particles  bouncing  off  of  the  road  surfaces  further  dis¬ 
turbing  the  dusty  road  substrate,  is  another  dust  generation  mechanism  that  probably 
plays  a  non- negligible  role.  Some  discussions  of  these  dust  generation  mechanisms  can 
be  found  in  the  literature  [26,  2],  The  amount  of  the  dust  emitted  into  the  air  by 
these  actions  depends  on  many  factors  such  as  dust  particle  properties  (size,  density, 
stickiness,  etc.),  environment  conditions  (wind  speed, etc.),  soil  conditions  (silt  content, 
moisture  level,  etc.)  and  vehicle  features  (vehicle  velocity,  weight,  dimension,  number  of 
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the  wheels,  etc.).  There  have  been  many  experimental  investigations  of  the  relationship 
between  the  magnitude  of  dust  emission  and  the  above  factors.  Experimental  meth¬ 
ods  include  utilizing  upwind/downwind  dust  flux  tower  measurement  system  [8,  7,  32], 
filter-based  gravimetric  techniques  [29,  30],  held  wind  tunnel  [18,  24,  6]  and  the  optical 
remote  sensing  method  [1].  The  U.S.  Environmental  Protection  Agency  (EPA)  has  col¬ 
lected  and  compiled  the  dust  emission  factor  information  gathered  under  different  held 
conditions.  The  data  and  dust  emission  models  are  documented  in  the  Compilation  of 
Air  Pollutant  Emission  Factors  (AP-42)  [29,  30,  28].  It  has  to  be  mentioned  that  most 
experimental  work  on  dust  emission  induced  by  vehicles  had  an  objective  to  model  dust 
emissions  as  line  sources  for  air  quality  modeling  and  so  near  held  collection  of  data 
was  not  a  priority.  For  our  application  in  which  the  vehicle  performance  is  the  main 
concern,  data  collected  close  to  the  vehicle  are  needed.  However,  experimental  data 
that  provide  detailed  near-vehicle  measurements  of  dust  concentrations  is  extremely 
limited  which  limits  the  extent  to  which  near  held  dust  emission  and  transport  models 
can  be  validated. 

As  for  modeling  of  dust  emission  and  transport  in  the  near-vehicle  region  there 
has  been  relatively  small  amount  of  work  described  in  the  open  literature.  Most  of 
this  work  focused  on  real-time  realistic  visualization  of  vehicles  such  as  the  numerical 
models  developed  at  George  mason  University [2,  3].  Their  model  focused  on  represent¬ 
ing  vehicles  in  a  virtual  environment  and  required  a  compuatationally  efficient  model 
since  interactive  graphics  required  the  model  operate  faster  than  physical  time.  As  a 
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result,  their  numerical  models  were  highly  simplified,  and  the  simulation  results  were 
not  validated  by  any  experimental  study. 

Turbulence  provides  another  challenge  in  the  modeling  of  dust  transport  around  the 
vehicle  once  it  is  emitted  into  the  air.  Once  in  the  air,  the  dust  is  scattered  by  small 
scale  turbulent  eddies  that  form  in  the  vehicle  wake  as  well  as  large  scale  turbulent 
eddies  that  are  superimposed  from  the  ever  present  atmospheric  boundary  layer  (ABL). 
Since  these  features  play  an  important  role  in  accurately  describing  the  mixing  and 
dispersion  of  dust  particles,  accurately  capturing  the  ABL  is  an  important  component 
of  modeling  dust  transport  about  moving  vehicles.  In  our  present  simulations  we  model 
the  turbulent  eddies  using  a  Reynolds  Averaged  Navier  Stokes,  RANS,  approach  and 
the  scattering  due  to  turbulence  is  modeled  as  a  stochastic  forcing  term  applied  to 
the  particles.  Care  must  be  taken  to  accurately  model  the  ABL.  In  particular,  it  can 
be  challenging  to  recover  an  ABL  profile  with  standard  engineering  RANS  turbulence 
models.  For  our  study  we  utilize  the  approach  that  has  been  shown  to  sustain  an  ABL 
profile  for  stable  atmospheric  conditions  over  a  long  fetch  [21,  33].  This  method  utilizes 
a  modified  surface  roughness  factor  along  with  an  imposed  stress  boundary  condition 
at  the  upper  boundary  to  create  a  correct  ABL  profile  that  is  sustained  over  long 
distances.  The  effects  of  turbulence  on  the  dispersion  of  dust  particles  have  been  well 
studied[12,  25,  19,  4,  22],  For  our  simulations  we  model  turbulent  fluctuations  using  a 
first  order  Markov  process  as  described  by  Zannetti[35]. 

We  utilize  the  verified  production  Loci/CHEM  multi-physics  CFD  code[31,  14]  to 


UNCLASSIFIED 


6 


simulate  the  dust  emission  and  transport  about  a  movin  vehicle  in  this  study.  The 
CHEM  code  is  a  robust  highly  parallel  solver  that  provides  a  variety  of  turbulence  mod¬ 
els  and  Lagrangian  particle  models  that  are  needed  to  model  the  physics  of  dusty  flows. 
In  addition  the  solver  provides  a  flexible  mechanism,  through  the  Loci  framework [36], 
for  extending  the  models  provided  by  the  baseline  CHEM  solver.  In  the  following  sec¬ 
tions  we  describe  the  development  of  specialized  models  for  the  emission  of  particles 
and  their  transport  and  dispersion  by  the  turbulence  in  the  vehicle  wake  and  ABL.  We 
study  the  sensitivity  of  the  model  to  source  locations  of  particles  as  well  as  particle  size 
distributions.  Finally  we  provide  a  careful  comparison  between  the  numerical  model 
and  experimental  data  providing  one  of  the  first  attempts  at  validating  a  numerical 
model  for  the  near  field  dust  transport  around  a  moving  vehicle. 

2  Numerical  Simulation  Models 

In  this  study,  numerical  simulations  were  carried  out  for  the  investigation  of  fugitive 
dust  emission  and  transport  around  a  moving  vehicle.  In  order  to  accurately  capture 
the  dispersion  of  the  fugitive  dust  emissions  the  turbulence  of  the  atmospheric  bound¬ 
ary  layer  must  be  considered.  When  using  engineering  level  RANS  based  turbulence 
models  special  care  must  be  taken  to  ensure  that  the  atmospheric  boundary  layer  is  sus¬ 
tained  over  the  simulation  domain.  The  failure  of  RANS  based  models  with  out  proper 
treatment  to  accurately  capture  the  ABL  has  been  described  in  the  literature [21,  33]. 
For  our  simulations  we  use  the  method  proposed  by  Richard  and  Hoxey  [20]  whereby 
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the  ABL  is  sustained  by  utilizing  an  appropriate  surface  roughness  parameter  on  the 
ground  along  with  an  imposed  stress  boundary  condition  on  the  upper  boundary.  We 
implement  these  schemes  using  Menter’s  baseline  (BSL)  turbulence  model [16]  and  in 
separate  studies  find  that  the  approach  is  effective  at  maintaining  the  ABL  over  the 
domain  sizes  considered  in  this  study  which  are  on  the  order  of  a  few  dozen  vehicle 
lengths. 

We  employ  the  Lagrangain  particle  tracking  method  to  compute  the  transport  of  the 
fugitive  dust  emissions.  Since  the  RANS  turbulence  models  do  not  resolve  the  turbulent 
eddies,  the  Lagrangian  particle  tracking  method  will  not  recover  the  dispersion  that 
these  eddies  cause.  In  order  to  correctly  reconstruct  the  scattering  effects  of  the  modeled 
turbulent  eddies  we  employ  a  stochastic  model  of  the  turbulent  fluctuations[34,  22, 
27].  In  particular  we  use  the  method  of  Zannetti[35]  whereby  the  turbulent  velocity 
fluctuations,  u',  are  modeled  as  a  simple  first  order  auto-correlation  Markov  process. 
In  this  formulation,  turbulent  fluctuation,  u',  is  expressed  as: 

u \t2)  =  R^h  -  tl)u'(ti)  +  u"(f2)  (1) 

where  Ri(t2—ti)  contains  the  auto-correlation  with  lag  At  =  t2—t\  of  the  u'  components 
where  u"  is  a  purely  random  vector.  LTsing  this  method  the  perturbation  of  each  particle 
is  independent  and  not  affected  by  the  position  of  other  particles  in  the  flow.  Therefore, 
the  numerical  performance  is  fast  since  no  special  interactions  between  particles  and 
fluid  are  required.  The  auto-correlation  of  the  Lagrangian  velocity,  RL  can  be  related 
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to  Lagrangian  time  scales,  Tl  by 

Rl  =  exp[-(t2  -  h)/TL\ 


(2) 


where  Lagrangian  time  scale  represents  the  time  over  which  the  velocity  of  a  particle  is 
self-correlated.  There  is  another  time  scale  -  Eulerian  time  scale  (Te)  that  represents 
the  time  the  turbulent  eddy  pass  through  a  fixed  point.  Normally  Lagrangian  time- 
scales  are  more  difficult  to  measure  as  it  requires  making  a  measurement  in  the  fluid 
frame  of  reference.  Many  models  have  been  suggested  to  ford  the  relation  between 
Lagrangian  and  Eulerian  time  scales  [23,  9,  10].  The  Lagrangian  time  scales  are  related 
to  the  Eulerian  time  scales  through  the  relation  Tl  —  r)TE  where  the  Eulerian  time 
scale,  Te,  can  be  derived  from  the  energy  dissipation  time  scale,  and  rj  is  assumed  to  be 
unity  from  other  experimental  results.  Assuming  that  turbulence  is  homogeneous  and 
isotropic  the  Eulerian  time  scale  can  be  determined  from  the  RANS  model  through  the 
specific  dissipation  rate  variable,  oj,  using  the  relation 


Te  = 


1  2 
3 


(3) 


where  C'/t  =  0.09.  Note,  if  we  were  to  consider  the  inhomogeneity  of  turbulence  in  lon¬ 
gitudinal,  lateral,  and  vertical  directions  then  the  1/3  factor  would  be  replaced  with  an 
appropriately  weighted  set  of  coefficients  for  each  of  the  corresponding  directions.  While 
the  assumptions  of  homogeneity  is  satisfactory  for  near  field  simulations  presented  here, 
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non- homogeneity  would  need  to  be  considered  for  more  distant  atmospheric  transport 
predictions  due  to  the  inhomogeneity  of  the  ABL. 

Next  term  to  be  modeled  in  equation(l)  is  the  random  vector  u".  It  is  assumed  that 
u"  follows  normal  distribution  of 


u"  = 


V2 


7TCT,," 


-.exp{- 


x 


2u2 


(4) 


where  cru«  is  the  standard  deviation  of  the  velocity  components,  characterized  by 


(Tn "  —  CU 


,\Ji-  -  RfjM  - 1\) 


(5) 


where  cru/  is  the  standard  deviation  of  u',  dehned  as 


where  k  is  the  turbulent  kinetic  energy. 


3  Modeling  Dust  Emissions 

For  the  purpose  of  simulating  the  near-held  generation  and  transport  of  dust  we  need 
models  for  the  quantity  of  dust  emitted,  the  location  of  the  emissions,  and  the  par¬ 
ticulate  size  distribution  of  the  emissions.  The  overall  mass  flux  of  particles  emit¬ 
ted  from  a  moving  vehicle  have  been  relatively  well  characterized  by  a  number  of 
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studies[30,  8,  7,  32,  13,  11,  17]  which  characterize  the  total  emissions  from  vehicle 
traffic  (suitable  for  use  in  line  source  models)  based  on  various  factors  such  as  vehicle 
speed  and  weight,  soil  silt  content,  and  soil  moisture  content.  These  studies  do  not  pro¬ 
vide  localized  information  about  the  distribution  of  dust,  but  studies  have  shown  that 
the  total  dust  concentrations  for  particles  smaller  than  10  microns  remain  relatively 
unchanged  over  distances  on  the  order  of  1  kilometer  provided  that  the  terrain  is  rela¬ 
tively  open[32].  Field  studies  show  that  the  quantity  of  dust  emission  depends  on  the 
parameters  that  characterize  the  condition  of  the  road  (e.g.  soil  silt  content,  moisture 
content),  vehicle  features  (e.g.  vehicle  weight,  vehicle  speed,  number  of  wheels)  and  the 
climate  conditions  (e.g.  frequency  and  amounts  of  precipitation).  Silt  content  refers 
to  the  mass  of  silt-size  material  (equal  to  or  less  than  75  /im  in  physical  diameter)  per 
unit  area  of  the  travel  surface.  The  downwind  dust  emission  measurements  in  emission 
models  were  usually  made  5  meters  from  the  road.  Emission  factor  (EF)  is  usually  the 
value  that  is  modeled  in  dust  emission  models,  and  it  is  usually  expressed  as  the  weight 
of  dust  divided  by  the  unit  weight,  volume,  distance,  or  duration  of  the  activity  that 
emits  the  dust.  For  air  quality  modeling  the  EPA  model  AP-42[28]  is  often  utilized 
to  model  particle  sources  clue  to  vehicle  traffic  on  unpaved  roads.  We  also  utilize  the 
AP-42  model  to  characterize  the  bulk  mass  flux  of  particle  emissions.  We  combine  this 
global  model  with  a  model  of  local  mass  emission  to  construct  a  near  field  dust  emission 
model.  The  components  of  this  model  are  outlined  in  the  following  subsections. 
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3.1  The  EPA  AP-42  Model 

We  adopt  the  emission  model  documented  in  AP-42  [28]  which  is  a  primary  compilation 
of  emission  factor  information  from  comprehensive  studies  conducted  by  the  U.  S.  Envi¬ 
ronmental  Protection  Agency  (EPA)  over  last  few  decades.  The  models  were  developed 
using  multiple  linear  regression  analysis  of  EF  test  data  versus  parameters  that  affect 
the  amount  of  particle  emissions.  Although  the  emission  factor  model  for  unpaved  road 
has  been  modified  over  the  past  20  years,  all  versions  have  important  common  features. 
The  silt  loading  has  consistently  been  found  to  be  the  most  important  factor  in  emission 
models.  In  this  study,  we  use  2003  version  of  AP-42  emission  factor  model  that  was 
developed  through  a  collection,  selection  and  analysis  of  data  from  12  field  test  reports 
on  dust  emission  induced  by  vehicles  traveling  on  unpaved  roads,  as  shown  in  eqn.  (7). 

^  Hs/unw/if 

EF  ~  (M/0.2)-  (7) 

where  EF  =  Emission  factor,  pounds  per  vehicle-mile-traveled,  (lb/VMT) 
k  =  Particle  size  multiplier  (dimensionless) 
s  =  Silt  content  of  road  surface  material  (%) 

W  =  Mean  vehicle  weight,  ton  (ton) 

M  =  surface  material  moisture  content  (%) 

For  air  quality  purposes  the  particulate  emissions  are  characterized  into  three  size 
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Table  1:  Parameters  for  eqn.  (7) 


constant 

PM2. 5 

PM10 

P  m30 

k 

0.38 

2.6 

10 

a 

0.8 

0.8 

0.8 

b 

0.4 

0.4 

0.5 

c 

0.3 

0.3 

0.4 

classifications  of  particle  matter  (PM)  known  as  PM30,  PM10,  and  PM2 5  which  rep¬ 
resent  all  particle  matter  that  has  a  diameter  less  than  30,  10,  and  2.5  micrometers 
respectively.  The  particle  size  multiplier,  k  for  different  size  range  for  eqn.  (7)  is  given 
in  table  1. 

Note  that  “normalizing  factors”  of  12  percent  silt  content  and  3  kilogram  vehicle 
weight  are  used  in  both  models.  0.2  percent  soil  moisture  content  is  selected  as  the 
reference  moisture  content  in  eqn.  (7),  which  is  usually  considered  as  the  default  dry 
condition  moisture  content.  The  moisture  content  is  influenced  by  meteorological  and 
physical  parameters  that  varies  with  time  and  location.  According  to  AP-42,  the  over¬ 
all  mean  moisture  content  in  publicly  accessible  road  data  set  was  1.1  percent.  The 
moisture  content  for  the  desert  ranges  from  0.17  to  0.48  percent.  Soil  silt  content  has 
significant  effect  on  dust  emission  factors.  A  summary  of  values  of  silt  content  on  in¬ 
dustrial  and  public  roads  is  given  in  table  13.2.2-1  in  AP-42  [28].  Mean  silt  contents  for 
gravel  roads  and  dirt  roads  were  found  to  be  6.4  percent  and  11  percent  respectively. 
Note  for  these  models  to  remain  consistent  with  experimental  conditions,  mean  wind 
speed  should  range  between  4  to  20  MPH. 
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3.2  Localized  Mass  Injection  Models 

The  EPA  AP-42  model  provides  a  reasonable  model  of  the  bulk  of  particulate  matter 
that  is  emitted  due  to  vehicle  travel  (at  uniform  speed  on  a  unpaved  road  surface), 
however  it  provides  no  information  about  where  the  emissions  emanate  from  a  vehicle. 
For  our  simulations  we  are  interested  in  the  detailed  transport  of  particulate  matter 
around  the  geometry  of  the  vehicle  and  for  this  we  need  to  know  the  actual  locations 
of  emissions.  Generally  emissions  will  come  from  three  sources:  1)  directly  from  the 
contact  points  (e.g.  tires  or  tracks)  of  the  vehicle  with  the  dusty  surface,  2)  saltation 
processes  caused  by  impacts  of  emitted  particles  with  the  surface,  and  3)  particles 
liberated  clue  to  the  action  of  shear  of  the  air  flow  created  by  the  vehicle  motion.  For 
dusty  unpaved  road  surfaces  we  assume  that  the  majority  of  dust  is  emitted  directly 
from  contact  points  rather  than  from  saltation  and  shear  processes,  which  appears  to 
be  the  case  from  simply  observing  vehicles  traveling  on  unpaved  roads. 

For  the  contact  point  emissions,  there  appears  to  be  two  primary  mechanisms  for 
emission:  1)  surface  adhesion  and  shear  with  the  ground  cause  particulates  to  stick 
to  the  tire  or  track  which  are  then  emitted  due  to  centrifugal  force  as  the  tire/track 
rotates  up  from  the  ground,  and  2)  crushing  action  at  the  leading  edge  of  the  tire/track 
causes  particles  to  be  ejected  at  the  leading  edge.  Unfortunately,  from  the  literature 
search,  we  have  not  found  any  experimental  data  that  provides  localized  information 
about  emissions  such  as  the  relative  quantity  of  particulate  matter  emitted  from  these 
localized  regions.  Lacking  quantitative  data,  we  have  created  a  model  that  qualitatively 
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produces  dust  plumes  near  the  tires  that  are  similar  to  observation. 

In  our  emission  model  we  assume  that  emissions  come  from  the  tire  tread  surface 
and  we  assign  a  probability  of  emission  as  a  distribution  over  the  tire.  In  this  model  we 
describe  the  tire  contact  with  the  ground  according  to  the  schematic  diagram  shown 
in  figure  1.  We  compute  an  emission  probability  based  on  the  angle  of  the  tire  tread 
surface  normal  with  respect  to  the  coordinate  system  formed  by  the  axle  vector  (into 
the  page)  and  the  gravitational  force  vector.  If  the  surface  normal  points  within  10 
degrees  of  the  axle  vector,  it  is  assumed  that  this  is  part  of  the  tire  sidewall  and  the 
emission  probability  is  set  to  zero.  Otherwise  the  probability  is  set  by  the  angle  6  as 
shown  in  figure  1.  The  non-normalized  emission  probability  is  given  by  the  equation 

£  =  e  *  (max (cos  6,  0})^  (8) 

where  e  is  an  efficiency  factor  that  is  used  to  scale  the  emissions  between  the  leading 
and  trailing  points  of  contact,  and  r/  is  a  parameter  that  describes  the  localization  of 
emissions.  The  maximum  function  prevents  emissions  from  occurring  on  the  top  of  the 
tire.  The  efficiencies,  e,  and  localization  parameter,  rj  for  the  trailing  contact  point 
(where  6  is  positive)  is  set  to  et  =  0.9  and  r}t  =  2,  while  the  leading  contact  point 
(where  6  is  negative)  is  set  to  ei  =  0.1  and  rft  =  6.  With  these  settings,  90%  of  the 
emissions  come  from  the  trailing  point  of  contact,  while  the  leading  point  of  contact 
has  a  distribution  of  emissions  that  is  heavily  biased  towards  the  ground  region.  The 
probability  distribution  function  is  formed  by  normalizing  £  by  the  integration  of  £ 
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Figure  1:  Schematic  of  Tire-Ground  Contact 

over  the  tire  surface.  Particles  are  assumed  to  be  emitted  from  this  surface  with  the 
tangential  velocity  of  the  rotating  tire.  We  use  this  model  to  describe  where  the  localized 
emissions  of  the  AP-42  model  emanate  from  in  the  vehicle  mesh,  while  the  overall  mass 
flux  is  adjusted  to  meet  the  AP-42  model  for  vehicle  dust  emissions. 

3.3  Qualitative  Evaluation 

Before  we  begin  validation,  we  perform  a  simulation  of  a  moving  vehicle  to  determine 
if  the  model  produces  qualitatively  reasonable  results.  For  this  test  case  we  consider 
the  simulation  of  a  Nissan  Pathfinder  SUV  for  which  geometry  was  already  available 
to  evaluate  the  simulation  of  dust  transport  for  a  realistic  configuration.  For  this  sim¬ 
ulation  we  utilize  the  Loci/CHEM  solver  using  a  15  million  cell  viscous  mixed  element 
mesh  that  as  generated  by  the  SolidMesh/AFLR3  mesh  generation  software  [15]. 

For  the  simulation  of  particle  generation  and  transport,  we  have  the  SUV  traveling 
with  13.5m/s  (about  30MPH)  uniform  forward  velocity  with  a  cross  wind  of  6.75  m/s 
(about  15MPH).  The  wind  is  simulated  as  an  ABL  with  a  velocity  of  6.75m/s  at  one 
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meter  height  using  a  surface  friction  velocity  of  w*  =  0.6 m/s  and  roughness  length  of 
y0  =  0.01m.  The  rubber  part  of  the  tire  is  prescribed  a  rotational  velocity  348.42  RPM 
which  is  consistent  with  a  forward  velocity  of  13.5m/s.  The  rotation  of  the  spokes  of 
the  hub  is  not  considered  in  this  model.  Particles  having  a  material  density  of  quartz 
(2648  kg/m3)  are  emitted  using  log-normal  particle  size  distribution  with  particle  mean 
diameters  of  30 ym  and  a  variance  of  a  =  0.25.  The  BSL  turbulence  model  and  the  first 
order  stochastic  turbulent  dispersion  model  are  employed  in  the  simulation.  Initial  dust 
mass  flux  released  by  tires  is  determined  using  model  described  by  eqn.(7).  The  silt 
content  and  soil  moisture  content  are  set  to  the  default  value,  namely  12%,  and  0.2% 
respectively.  The  weight  of  vehicle  is  taken  from  Nissan  Pathfinder,  which  is  4500 kg. 
The  mass  flux  computed  from  the  model  is  then  equally  distributed  to  the  4  wheels. 
Particles  are  ejected  from  the  tires  using  the  localized  mass  injection  model  represented 
by  eqn.(8). 

Side  and  top  views  of  the  simulated  dust  distribution  around  the  vehicle  for  this 
case  are  shown  in  figure  2.  The  vehicle  surface  and  ground  are  shaded  by  gage  pressure 
contours.  Instead  of  viewing  particles  as  points,  we  use  the  volume  visualization  method 
in  which  the  computed  particle  concentration  field  is  averaged  over  the  volume  in  time 
allowing  the  particle  cloud  to  be  viewed  as  a  semi-transparent  region  to  mimic  the 
scattering  effects  that  visually  dominate  particle  plumes.  From  this  figure  it  can  be 
seen  that  the  localized  release  model  of  particles  produces  a  distribution  of  dust  both 
around  the  tires  and  the  vehicle  that  follows  qualitative  expectations  (e.g.  the  simulated 
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particle  cloud  looks  similar  to  what  is  observed  when  a  vehicle  drives  down  a  dusty  road). 
Particles  are  initially  released  near  the  tire  contact  points  and  are  then  entrained  by  the 
recirculating  flows  in  the  vehicle  wake  and  scattered  by  both  the  wake  and  atmospheric 
turbulence.  The  structure  of  the  mean  wake  of  the  vehicle  is  not  strongly  evident 
because  turbulent  dispersion  plays  a  role  on  scattering  the  particles.  Note,  that  since 
the  stochastic  model  does  not  include  any  spatial  correlations,  there  is  no  billowing 
type  of  effect.  Each  particle  follows  its  own  “Brownian  Motion”  type  of  trajectory,  but 
a  group  of  particles  in  the  same  region  of  space  do  not  follow  the  same  “eddy”  and  the 
overall  effects  are  averaged  out  to  give  an  average  characterization  of  the  particle  plume. 
In  general,  these  results  show  that  the  described  model  performs  in  a  fashion  that  is 
qualitatively  comparable  to  the  observed  dusty  cloud  produced  by  vehicles  traveling  on 
unpaved  surfaces. 
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Figure  2:  Dust  distribution  around  the  moving  vehicle  from  views  on  front  driver’s  side 
(top),  front  passenger  side  (middle)  and  top  (bottom)  with  volume  visualization 
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4  Model  Validation  Case 

In  order  to  validate  the  model  we  would  like  to  compare  to  experimental  data  that  is 
as  close  to  the  vehicle  as  possible  in  order  to  better  validate  near  held  particle  distri¬ 
butions.  Unfortunately,  most  experimental  data  are  collected  at  some  distance  from 
the  vehicle  since  their  purpose  is  to  validate  environmental  quality  models  that  treat 
dusty  roads  as  line  sources.  After  reviewing  the  available  literature,  we  have  deter¬ 
mined  that  measurements  of  dust  concentrations  carried  out  by  Veranth  et.al  [32,  5] 
provide  the  most  appropriate  data  for  our  validation  needs.  The  main  reason  we  chose 
this  data  set  is  that  there  are  measurement  data  conducted  at  3  meters  away  from  the 
road,  while  other  measurements  from  our  literature  survey  were  performed  at  much 
greater  distances.  In  addition,  the  data  presented  in  these  experiments  included  ab¬ 
solute  measurements  of  dust  concentrations  that  are  suitable  for  validation  exercises, 
instead  of  relative  measurements  that  were  designed  to  calibrate  existing  models.  Fi¬ 
nally,  compared  to  other  experimental  work,  the  running  conditions  in  this  held  study 
were  relatively  well  documented  which  improves  the  suitability  of  this  data  set  for 
validating  our  detailed  near-held  models. 

4.1  Outline  of  Experimental  Data  Set 

The  held  measurements  were  conducted  at  the  Dugway  Proving  Ground,  Tooele  County, 
Utah  by  University  of  Utah  in  collaboration  with  the  Mock  Urban  Setting  Test(MUST). 
The  original  purpose  of  this  study  was  to  investigate  the  signihcance  of  near-source 
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removal  of  the  dust  under  the  stable  atmosphere  and  high  surface  roughness  conditions. 
Surface  roughness  that  might  be  presented  by  large  surface  roughness  like  buildings  was 
modeled  in  the  experiment  by  a  10  by  12  array  of  2.5m  high,  2.4m  wide,  and  12.2  long 
rectangular  cargo  shipping  containers.  In  the  surrounding  landscape,  vegetation  is  a 
thin  cover  of  brush  of  0.5-lm  high.  The  soil  was  classified  by  the  Natural  Resources 
Conservation  (2000)  as  Skumpah  silt  loam  and  had  16%  silt  content.  The  atmospheric 
condition  was  stable.  The  test  was  performed  using  a  1994  Ford  pickup  truck  with 
3900 kg  weight  that  was  driven  at  the  speed  of  9 m/s  on  an  180m  length  of  unpaved  road 
that  ran  parallel  to  the  upwind  edge  of  the  container  array.  Wind  speed  was  normally 
between  1  and  5 m/s  and  the  wind  direction  remained  within  45  degree  perpendicular 
to  the  road.  Table  2  shows  the  measured  wind  speed  at  different  height  with  standard 
deviation  averaged  within  1.5  hour  test  period.  Average  air  temperature  and  pressure 
were  294 K  and  0.855 atm  respectively. 

Table  2:  Mean  and  standard  deviation  wind  speed  for  1.5  hour  test  period 


Height  (m) 

Upwind  wind  speed  (m/s) 

4 

3.42  ±0.51 

8 

4.07  ±0.50 

16 

4.72  ±  0.52 

Dust  concentration  was  measured  using  DustTrak  analyzers  with  PM10  inlets.  The 
instruments  were  mounted  on  a  movable  tower  at  various  heights  to  collect  the  data  at 
different  location.  As  we  are  interested  in  the  near-field  predictions,  we  focus  only  on 
the  near-source  tower  location  which  is  located  3  meters  from  the  edge  of  the  road,  or 
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4.5  meters  from  the  road  center-line.  Dust  measurements  were  performed  on  the  tower 
at  heights  of  0.9,  1.7,  and  3.7  meters  above  the  ground.  The  vehicle  was  traveling  for 
total  1.5  hours  as  a  series  of  44  trips  at  1-1.5  minute  intervals,  and  the  time-averaged 
dust  concentration  was  recorded  every  5  seconds.  Peak  concentration  value  for  each  trip 
(mg/m3)  and  pulse  area  defined  as  the  time  integrated  concentration  per  trip  (mgs/m3) 
were  recorded  in  the  experiment  and  are  provided  in  table  3. 


Table  3:  Mean  and  standard  deviation  of  dust  concentration  data  at  near  source  (3 
meters  from  the  road) 


Height  (m) 

peak  concentration  (mg/m3) 

pulse  area  (mgs/m3)  per  trip 

0.9 

38.9  ±21 

302  ±  171 

1.7 

19.9  ±  11 

144  ±  92 

3.7 

10.3  ±  7.2 

77.8  ±  56 

4.2  Simulation  Setup 

For  the  simulation  setup  we  perform  simulations  in  the  vehicles  reference  frame.  In  this 
reference  frame,  the  tower  moves  through  the  domain  at  the  rate  of  the  vehicle  speed, 
and  we  can  obtain  the  time  history  of  the  tower  particle  concentrations  by  plotting  the 
concentration  along  a  line  that  represents  the  measurements  at  a  given  tower  height. 
In  this  plot,  distance  along  the  line  can  be  converted  into  time  in  the  tower’s  reference 
frame  simply  by  dividing  by  the  vehicle  speed.  Thus  we  have  a  straightforward  way 
to  compare  simulations  to  measured  experimental  data  provided  that  the  simulation 
domain  is  large  enough  that  the  entire  pulse  of  the  particle  plume  is  recorded  within 
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the  domain.  The  first  step  in  this  process  is  the  generation  of  a  mesh  for  a  long 
domain  that  included  the  vehicle  and  enough  downwind  region  to  properly  simulate 
the  plume’s  intersection  with  the  measurement  tower.  In  the  experiment,  a  1994  Ford 
pickup  truck  was  used  as  the  dust  emission  source.  For  the  geometry  of  the  1994  Ford 
pickup  we  constructed  a  simplified  model  based  on  general  vehicle  geometry  available 
in  generic  descriptions  from  Ford  including  vehicle  wheel-base,  bed  dimensions,  and  cab 
height.  Figure  3  shows  the  resulting  simplified  pickup  truck  geometry.  To  determine 
the  sensitivity  of  the  dust  plume  at  the  first  tower  location  to  the  vehicle  geometry  we 
ran  simulations  using  the  detailed  Pathfinder  SUV  and  the  approximate  pickup  truck 
geometry  and  found  that  there  were  not  significant  differences  in  the  simulated  dust 
plume  at  the  3  meter  tower  location.  This  is  not  surprising  as  it  confirms  observations 
in  the  work  of  Gillies  et.al  [8]  that  showed  that  the  size  of  the  wake  created  by  a  moving 
vehicle  was  dependent  on  the  height  of  the  vehicle.  In  our  study,  the  heights  of  SUV 
and  pickup  truck  is  about  the  same,  namely  1.7m  and  1.9m  respectively.  From  these 
studies  we  conclude  that  the  simplified  geometry  of  the  pickup  truck  is  sufficient  for 
validating  to  the  tower  data  that  is  3  meters  from  the  edge  of  the  road. 

The  experimental  data  setup  that  we  are  using  for  validation  included  an  array  of 
cargo  containers  that  were  placed  about  6  meters  from  the  road  with  the  first  dust  tower 
being  just  3  meters  away  from  the  road  edge.  Although  the  dust  tower  is  upwind  of  the 
container  array,  the  effect  of  the  blockage  represented  by  the  container  array  will  cause 
an  uplifting  flow  caused  by  the  blockage  of  the  array.  To  correctly  capture  the  effects  of 
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Figure  3:  Geometry  and  surface  mesh  of  a  Ford  Pickup  truck 

the  container  array  in  our  validation,  we  create  a  rectangular  block  with  the  same  height 
of  the  containers  (2.5 m)  located  at  6 m  from  edge  of  the  road  (distance  is  estimated  from 
the  schematics  given  in  the  report  [32])  and  extending  along  the  road  throughout  the 
domain  as  shown  in  figure  4.  The  measurement  tower  is  placed  half  way  between  the 
blockage  and  the  edge  of  the  road.  Simulations  that  we  have  performed  comparing  dust 
plumes  with  and  without  this  blockage  showed  that  it  contributes  a  significant  effect 
in  the  height  of  the  dust  plume  at  the  closest  measurement  tower.  In  our  simulation 
we  approximate  the  container  array  as  if  it  is  a  continuous  obstruction  along  the  road 
which  approximates  the  blockage  effect  of  the  containers  while  also  facilitating  a  simple 
steady-state  simulation  procedure. 

For  particle  surface  interactions  we  assume  that  any  particle  that  interacts  with  the 
vehicle  surfaces  or  with  the  ground  rebounds  specularly  off  of  the  surface.  Although 
some  particles  will  be  removed  on  contact  with  the  ground,  empirical  evidence  suggests 
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Figure  4:  Computational  domain  with  a  blockage  parallel  to  the  vehicle  traveling  di¬ 
rection.  Yellow  part  is  the  blockage,  and  particles  are  represented  by  red  color 

that  in  the  near  field  only  a  very  small  amount  of  particles  will  be  removed  from  the 
simulation  by  way  of  deposition.  Therefore,  we  simply  assume  that  no  particle  depo¬ 
sition  occurs  on  these  surfaces.  For  simulations  over  larger  distances,  a  proper  ground 
deposition  model  may  be  required  to  properly  consider  ground  particle  interactions. 

Dust  is  emitted  from  the  tires  with  a  distribution  that  is  described  by  equation 
(8)  with  the  mass  emission  equally  distributed  to  the  four  tires.  The  total  mass  flux 
of  the  tire  emissions  are  determined  by  the  AP-42  model  as  described  by  equation(7). 
The  parameters  of  16%  silt  content,  3900kg  vehicle  weight,  0.2%  moisture  content  with 
particle  size  multiplier  from  PM10  were  input  to  the  eqn.(7)  to  get  dust  emission  factor. 
This  value  combined  with  the  vehicle  speed  of  9 m/s  gives  0.0096 kg/s  total  mass  flux 
injected  from  each  of  the  four  tires.  A  model  for  the  distribution  of  the  particle  sizes  is 
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prescribed  by  a  log-normal  distribution  with  a  mean  particle  size  of  10/im  with  standard 
deviation  of  0.35. 

To  simulate  the  atmospheric  wind  effects,  a  background  atmospheric  wind  profile  is 
prescribed  using  logarithmic  law  as  follows: 


u*  ,z  +  z0^ 

u  =  —ml - 

A  zO  ' 


(9) 


where  K  =  0.4,  and  friction  velocity  ( u *)  and  roughness  height  (zo)  are  estimated  to  fit 
the  upwind  data  within  the  variation  of  measurements  listed  in  table  2.  u*  =  0.24 m/s 
and  Zo  =  0.01m  are  selected  in  the  wind  profile  expression.  The  value  of  friction  velocity 
agrees  with  held  measurement  of  from  sonic  anemometer  data  ( u *  =  0.23m/s)  very 
well.  Figure  5  demonstrates  that  the  prescribed  wind  speed  profile  compares  well  with 
experimental  data.  As  mentioned  in  held  study,  wind  direction  was  within  45  degrees 
perpendicular  to  the  road,  but  there  was  no  data  on  average  wind  direction.  In  our 
simulation,  the  wind  direction  was  hrst  set  perpendicular  to  the  road,  then  45  degrees 
to  the  road  to  evaluate  the  bounding  effects  of  wind  direction  in  the  validation. 


4.3  Validation  Results 

Simulations  for  two  different  wind  directions  (90°  and  45°)  were  conducted  for  com¬ 
parison  to  experimental  data  at  the  closest  tower  location.  The  dust  concentration 
temporal  history  for  each  of  the  three  tower  sampling  heights  (0.9,  1.7,  and  3.7  meters) 
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Figure  5:  Comparison  of  prescribed  atmospheric  boundary  layer  profile  and  experimen¬ 
tal  data  with  variation  of  measurement 
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were  determined  from  the  steady  state  simulation.  These  tower  heights  are  represented 
in  the  figure  6  showing  three  lines  that  represents  the  three  sampling  heights  drawn 
on  a  cutting  plane  located  at  the  tower  location.  From  these  spatial  lines  a  temporal 
history  can  be  reconstructed  as  shown  in  figure  7.  Integrating  the  temporal  histories 
we  determine  the  peak  and  integrated  concentrations  of  particles  at  each  of  the  tower 
heights  for  comparison  with  experimental  data  as  is  shown  in  table  4. 

The  simulated  peak  dust  concentration  compare  reasonably  well  with  experimental 
data  and  are  generally  with  experimental  error  bounds.  We  note  that  the  wind  direction 
used  in  the  simulation  has  a  significant  impact  on  the  resulting  dust  distribution  at  the 
tower  location.  In  general,  the  simulations  at  the  larger  wind  angle  appear  to  be  more 
consistent  with  the  experimental  data,  and  the  simulations  appear  to  under-predict 
integrated  dust  concentrations  suggesting  that  perhaps  more  particles  should  have  been 
emitted  from  the  vehicle  tires.  In  addition,  it  appears  that  the  experimental  results 
showed  a  slightly  more  dispersed  particle  plume  than  the  simulation,  in  particular  when 
observing  the  comparison  of  the  integrated  particle  concentrations  at  the  highest  tower 
sampling  location.  However,  our  simulations  are  still  within  the  experimental  error 
bounds  and  so  no  definitive  conclusion  can  be  made.  From  this  simulation,  it  is  also  clear 
that  the  wind  characterization  is  a  significant  source  of  variability  in  the  measured  and 
simulated  data,  thus  contributes  uncertainty  in  the  process  of  validation/calibration. 

As  already  mentioned,  numerical  studies  showed  that  details  of  the  vehicle  geometry 
did  not  play  a  significant  factor  in  the  measured  dust  plume.  In  addition  we  performed 
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Table  4:  Comparison  of  dust  concentration  results  between  CFD  simulation  and  ex¬ 
periment  at  closest  tower  (3m  from  road)  with  perpendicular  and  45°  simulated  cross 
winds _ 


Height 

(m) 

measured 

peak  concentration 
(mg /m3) 

simulated 

peak  concentration 

(mg/m3) 

measured 
pulse  area 
(mg  s/m3) 

simulated 
pulse  area 
(■ mg  s/m3) 

90° 

45° 

O 

O 

45° 

0.9 

38.9  ±21 

52.2 

31.0 

302  ±  171 

116. 

98.8 

1.7 

19.9  ±  11 

48.0 

26.0 

144  ±  92 

94.0 

84.6 

3.7 

10.3  ±7.2 

10.1 

7.51 

77.8  ±  56 

17.0 

27.6 

a  number  of  additional  simulations  to  determine  the  effect  of  other  factors  that  may 
effect  comparisons  to  the  numerical  experiments.  The  first  consideration  is  the  effects 
of  the  size  distribution  of  the  emitted  dust  particles.  For  the  dust  size  distribution  the 
amount  of  dust  in  the  different  categories  can  be  inferred  from  the  EPA  model  which 
provides  different  coefficients  for  PM2. 5,  PM10,  and  PM30  dust  size  categories.  From 
this  we  inferred  both  the  corresponding  log- normal  and  Weibull  size  distribution  PDFs 
and  compared  simulations  using  either  log-normal  (as  used  in  the  presented  results) 
and  the  Weibull  distribution  which  had  a  stronger  bias  towards  small  particle  sizes.  We 
found  that  there  was  no  significant  difference  in  the  measured  dust  plume  between  these 
two  PDF’s  suggesting  that  particle  sizes  played  a  minor  role  in  the  plume  dispersion  at 
the  nearest  tower.  In  addition,  we  evaluated  the  effects  of  the  exponent  on  the  localized 
dust  emission  model  given  in  eqn.  (8)  by  comparing  a  linear  distribution  (77  =  1)  instead 
of  the  square  law  used  in  the  validation  study.  In  terms  of  the  plume  shape  observed  at 
the  3  meter  tower  location,  changes  in  the  tire  release  model  did  not  significantly  change 
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Figure  6:  Dust  distribution  on  the  near  source  plane  (3m  from  the  road) 


the  dust  distribution.  However,  the  exponents  that  we  used  in  the  study  produced  a 
dust  plume  that  better  matched  the  qualitative  character  of  tire  emissions  observed 
from  vehicles  moving  on  dusty  roads. 
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Figure  7:  The  simulated  dust  concentration  pulse  for  different  wind  direction  at  0.9m 
(top),  1.7m  (middle)  and  3.7m  (bottom)  measuring  stations. 
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5  Concluding  Remarks 

Many  aspects  of  vehicle  functionality  depend  on  the  emission  and  transport  of  fugitive 
dust  caused  by  vehicle  movement,  especially  when  traveling  on  unpaved  surfaces.  Given 
the  needs  for  particle  emission  models,  there  are  very  few  detailed  models  of  the  fugitive 
dust  in  the  vehicle  near-held  region.  In  this  paper  we  presented  a  modeling  approach  for 
fugitive  dust  emissions  that  can  capture  near-held  dust  concentrations.  We  validated 
this  model  by  comparing  simulated  plumes  at  3  meter  towers  with  experimental  data 
and  found  that  the  proposed  model  performed  within  experimental  error  bounds.  This 
work  also  provided  the  hrst  attempt  to  validate  the  numerical  models  for  the  near  held 
dust  emission  and  transport  around  a  moving  vehicle. 

One  unambiguous  hireling  of  literature  review  and  our  simulation  studies  is  that 
current  publicly  available  experimental  data  is  inadequate  to  fully  validate  near-held 
dust  emission  and  transport  models.  Most  experimental  data  was  commissioned  to 
support  atmospheric  dust  transport  modeling  for  environmental  quality  impact  studies. 
Basically  these  experimental  studies  were  targeting  data  that  is  assumed  to  be  averaged 
over  time  periods  in  the  range  of  weeks  to  months,  and  over  length  scales  on  the  order 
of  kilometers  or  greater.  As  a  result,  most  of  these  studies  did  not  provide  definitive 
data  on  the  dust  in  the  local  area  of  the  vehicle  in  the  time-scale  of  a  single  traversal  of 
a  dusty  terrain.  Lacking  detailed  experimental  data  on  the  vehicle  spatial  and  temporal 
scales,  the  ability  to  refine  models  beyond  what  is  presented  in  this  study  is  somewhat 
limited:  while  more  sophisticated  modeling  approaches  could  be  suggested,  there  would 
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be  no  way  to  provide  sufficient  validation  of  these  models  based  on  data  available  in  the 
literature.  Therefore,  further  enhancement  of  these  models  must  proceed  along  with 
detailed  experimental  work  focused  on  near  field  measurement  of  fugative  dust  emission 
on  the  appropriate  time  and  space  scales  that  would  accurately  capture  localized  dust 
behavior. 
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